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Abstract 

We study muon-spin rotation (/iSR) spectra in the mixed phase of highly 
anisotropic layered superconductors, specifically Bi2+xSr2-xCaCu20g,^s 
(BSCCO), by modeling the fluid and solid phases of pancake vortices using 
liquid-state and density functional methods. The role of thermal fluctua- 
tions in causing motional narrowing of ^SR lineshapes is quantified in terms 
of a first-principles theory of the flux-lattice melting transition. The effects 
of random point pinning are investigated using a replica treatment of liq- 
uid state correlations and a replicated density functional theory. Our results 
indicate that motional narrowing in the pure system, although substantial, 
cannot account for the remarkably small linewidths obtained experimentally 
at relatively high fields and low temperatures. We find that satisfactory 
agreement with the /iSR data for BSCCO in this regime can be obtained 
through the ansatz that this "phase" is characterized by frozen short-range 
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positional correlations reflecting the structure of the liquid just above the 
melting transition. This proposal is consistent with recent suggestions of a 
"pinned liquid" or "glassy" state of pancake vortices in the presence of pin- 
ning disorder. Our results for the high-temperature liquid phase indicate that 
measurable linewidths may be obtained in this phase as a consequence of den- 
sity inhomogeneities induced by the pinning disorder. The results presented 
here comprise a unified, first-principles theoretical treatment of /iSR spectra 
in highly anisotropic layered superconductors in terms of a controlled set of 

approximations . 

PACS: 74.60.-w, 74.72.Hs, 76.75.+i 
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I. INTRODUCTION 

In the mixed phase of a type-II superconductor, an externally applied magnetic field 
penetrates the bulk of the sample in the form of lines of magnetic fiux [|l| . The distribution 
function n{B) of the local magnetic induction B in the sample is thus a non-trivial quan- 
tity, depending on the fiux distribution associated with a single vortex line, as well as on 
the arrangement of the vortex lines in space and time. This distribution function can be 
obtained through muon-spin rotation (/iSR) experiments based on the following procedure: 
An ensemble of muons, initially polarized transverse to the magnetic induction, is implanted 
uniformly throughout the sample 1|,|,H . Since the local magnetic induction varies in space, 
the spins of the muons in different regions of the sample precess at different rates, resulting 
in a dephasing of the initial polarization with time. This dephasing is monitored by track- 
ing the angular distribution of the positrons emitted as the muons decay. The measured 
polarization can then be related, via a Fourier transform, to n{B). 

The utility of the /iSR method is that it provides an accurate local probe of magnetic 
field variations in the bulk irrespective of whether such variations are periodic (as for a 
regular vortex lattice) or have more complex structure refiecting arrangements of vortex 
lines with varying degrees of spatial correlations. In high-temperature superconductors such 
as Bi2+xSr2-xCaCu208+s (BSCCO), /iSR experiments carried out in a field perpendicular 
to the superconducting layers (this is the geometry we consider throughout this paper) have 
provided indirect evidence for what is believed to be a melting transition |^ of the Abrikosov 
fiux-line lattice to a fiux liquid with short-ranged positional correlations. In the experiments 
of Lee et al [P,^^ on BSCCO, where n{B) is obtained as a function of the external field 
H and the temperature T, the location of the melting line at fields of the order of 10 mT 
is inferred from the abrupt changes in the /iSR linewidths and lineshapes which occur as 
H and/or T are varied across the transition boundary. In this low-field regime, the results 
obtained are in broad agreement with extensive experimental and theoretical work over the 
past decade aimed at detecting and understanding the nature of such a melting transition. 



In contrast, the nature of the transition and of the low-temperature phase at high fields still 
remains controversial. 

For a regular Abrikosov lattice at zero temperature and in the absence of disorder, 
the second moment of n{B) can be simply related to the magnetic penetration depth A 
1^. For a wide variety of superconductors, the values of A obtained in this way P, pil|JTT[| 



agree closely with those obtained by other means ||I2|. However, for BSCCO, such an 
interpretation of the high field data yields values for the ab plane penetration depth Xab 
ranging between 3500 and 4500 A |jl3|. These numbers are to be contrasted to typical 



values obtained by other means which suggest 1100 < Xab < 2800A, with values in the range 
1500 < Xab < 2000A believed to be appropriate for slightly overdoped samples R||,0,|T5|,|16| . 



This discrepancy originates in the anomalously small values of the linewidth J[AB'^] (the 
brackets [■ ■ ■] denote a space average) obtained in experiments performed at high fields [ll7| . 
Fits to the lineshapes at low fields, in contrast, yield results in reasonable agreement with 
conventionally accepted values for Xab with Xab — I8OOA obtained in the recent experiments 
of Lee et al. [^,^,§1 . The magnitude of the linewidths is found to change rather sharply across 
a "critical" value Her (~ 100 mT) of the magnetic field which depends only weakly on the 
temperature (the "low" and "high" field regimes referred to earlier can be more precisely 
correlated to Her)- As the external field is increased across Her-, the /iSR lineshapes are 
observed to become far more symmetric, concomitant with the sharp reduction in linewidths. 
This asymmetry can be quantified through the normalized "asymmetry parameter" a = 
[AB^Y^'^/IAB'^Y^'^. In the high field regime, measured values of a range from 0.5 - 0.9 
{a ~ 1.2 in the Abrikosov flux lattice phase). These observations suggest that the nature 
of the flux array changes qualitatively as the applied field is increased across Her- The 
existence of such a critical field value across which the attributes of the fiux-line system 
change abruptly has been confirmed in a variety of other measurements, including transport 
||T8| and magnetization [|I^ experiments and neutron scattering |^. Possible explanations 
for these observations include the proposal of a sharp crossover from three-dimensional to 



effectively two-dimensional behavior as a function of the field in the absence of disorder ^1 



an entanglement transition of the flux-line lattice [0,0 as the fleld is increased, and the 



transition between a pinning-induced, weakly disordered "Bragg glass" phase [0 at low field 
values 12^J2^ to a strongly disordered "pinned liquid" or "vortex glass" [0 phase at high 
fields. Figs. |V11| -P illustrate schematic phase diagrams. The phase diagram of a pure system 
according to the conventional Abrikosov theory is shown in Fig. |V11[ . A phase diagram 
illustrating the Abrikosov lattice and vortex liquid phases expected for a pure system is 
shown in Fig. |^, and Fig. ^ shows a conjectured phase diagram in the presence of random 
point pinning, illustrating Bragg glass, pinned liquid (or vortex glass) and vortex liquid 
phases. 

The explanations proposed for reduced linewidths in the high field regime of BSCCO fall 
into two major classes which are not completely unrelated. One class of theories suggests 
that reduced linewidths arise due to thermal broadening (and possibly quantum broadening 
- we ignore quantum effects in the bulk of this paper, but discuss them briefly in the 
concluding section) of density distributions in the lattice phase. It is argued that the effects 
of thermal fluctuations in BSCCO are enhanced above H^r due to the nearly two-dimensional 
character of vortex fluctuations in this regime - the flux-line structure in BSCCO is more 



accurately thought of in terms of "pancake" vortices [^ moving on the superconducting 
layers and interacting via a combination of electromagnetic and Josephson couplings. Such 
behavior should generalize to all sufficiently anisotropic, layered compounds ||2^ and should 
be refiected to a lesser extent in more isotropic materials. Calculations representing this class 
of theories have so far been restricted to harmonic treatments of elastic fiuctuations about 
the zero-temperature Abrikosov fiux lattice 0J|,0,0]. In these calculations, the effects of 
thermal fluctuations in the lattice phase are expressed in terms of the Lindemann parameter 



L = J{u'^)/aQ, where (u^) is the mean square thermal fluctuation of a flux-line from its 
equilibrium position at melting and Qq is the mean inter-line spacing. In the analysis of 
experimental data, L is treated as a phenomenological fltting parameter. One purpose of 
this paper is to present calculations of /iSR lineshapes and linewidths based on a detailed 
and microscopic theory of flux-liquid and flux-lattice phases in a highly anisotropic layered 



superconductor. This theory can be used to calculate the Lindeniann parameter as well as 
density distributions at freezing. Such density distributions feed into a calculation of the 
lineshapes measurable in /xSR experiments - results of a calculation of these lineshapes are 
presented here. 

The second class of theories invokes the presence of underlying quenched disorder and pro- 
poses random pinning as a possible source of reduced linewidths. The quasi-two-dimensional 
vortex system above Her is argued to be highly susceptible to pinning. Several proposals 
for the nature of this low-temperature state exist in the literature. The "vortex glass" de- 



scription [^ concentrates on the physics at large length scales and thus does not address 



issues of short-range order. There has been much recent interest in the proposal ||2^ that 



at low fields, the flux-line system is in a "Bragg glass" phase in which the (5-function Bragg 
peaks associated with the pure crystalline phase are replaced by algebraic singularities. 



Such a relatively ordered phase is expected p^j25|j26[] to yield to a topologically disordered 
phase (where dislocations proliferate) at fields larger than a critical value, this transition 
being driven by the increasing relevance of disorder due to the effective reduction of the 
dimensionality at high fields. Evidence for such a transition has been found in numerical 



simulations |^ of a system of flux lines in the presence of random point pinning. The 
Bragg glass scenario is a plausible explanation of the experimental data. However, the pre- 
cise nature of the phase obtained at large fields remains unclear. Given that the physics 
of the problem at high fields is dominated by the fluctuations of pancake vortices weakly 
coupled across layers, correlations of vortices across the superconducting layers are expected 
to be much weaker than correlations between vortices on the same layer. Specifically in 
the context of /iSR experiments on BSCCO, some analytic work exists on a model which 
assumes that the low-temperature state in the presence of pinning disorder is characterized 



by a disordered stacking of ordered planes p0|,|32| . However, it is difficult to see how even 



relatively strong disorder could stabilize such an arrangement in which the misalignment of 
ordered layers is penalized by a macroscopic (proportional to the area of the layer for each 
layer) energy cost. It is reasonable to expect that some fraction of this cost is relaxed by 
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allowing dislocations to form on each layer. If a non-zero density of dislocations is allowed 
to form, the detailed microscopic arrangement of vortices must be qualitatively different 
from that proposed in such simplified models. In this paper we propose that a reasonable 
assumption for this disordered state is that it resembles a "frozen" liquid in the sense that 
the positional correlations in this state are essentially the same as that in a liquid just above 
its putative freezing transition, but the time scale for structural relaxation is much longer 
than that in the liquid. 

In our theoretical analysis of /iSR spectra in the mixed phase of extremely anisotropic, 
layered, type-II superconductors in a field applied perpendicular to the superconducting 
layers, the flux-solid phase is described using the density functional theory |]33| , |3^ , |35| , |36| and 
the fluid phase using a liquid state-theory adapted by us p5|j36|l to study "pancake" vortices 
in the Lawrence-Doniach description [^ of a layered superconductor in the limit of vanishing 
Josephson couphng. Our approach assumes a value of \ab consistent with the values obtained 
from low-field /iSR and other experiments on BSCCO (1500A < Xab < I8OOA) and calculates 
linewidths and lineshapes appropriate to these values, incorporating thermal fluctuations, 
disorder effects (we consider random point pinning produced by atomic-scale pinning centers 
such as oxygen defects [^) as well as corrections due to the finite core size [^. The 
approach outlined here should be applicable to the calculation of /iSR lineshapes in artificial 
superconducting multilayers [Q in which the distance between neighboring superconducting 
layers can be tuned so as to reduce the Josephson coupling between layers to an insignificant 
level. The approach of this paper is, in principle, non-perturbative in that our description of 
thermally broadened density distributions in the flux-solid phase is not restricted to assuming 
a harmonic form for the free energy cost associated with small deviations from the perfect 
lattice. 

Our principal results for the pure system include a calculation of yuSR lineshapes and 
linewidths at freezing. These results can be extended below the freezing transition by as- 
suming harmonic fluctuations at sufficiently low-temperatures, whose magnitude is consis- 
tent with the theoretically calculated (m^) at the melting transition. We find that thermal 



broadening should reduce the /iSR hnewidth for BSCCO at the flux-lattice melting temper- 
ature by a factor of about 4 relative to that for a perfect Abrikosov lattice. These results 
quantify thermal fluctuation amplitudes in the pure system and the magnitude of the cor- 
responding corrections to the linewidths. These corrections can then be used to estimate 
Xab more accurately. However, the linewidths so obtained are still far in excess of those 
measured experimentally for H ^ Her even at the "melting" transition. Further, theories 
based on thermal fluctuations about an ordered triangular lattice of flux lines cannot explain 
the anomalously small linewidths obtained at low temperatures, as pointed out by several 
other authors. Thus, disorder effects appear to be necessary for an explanation of the /xSR 
linewidths at high fields. 

We also study, using a replicated liquid-state and density functional theory [^, the 



effects of weak pinning disorder on the /iSR linewidth in both the high-temperature liquid 
and the low-temperature solid. We find that such disorder should not affect /iSR lineshapes 
considerably in the low-temperature state if it has strong local crystalline order, such as 
might be expected in a state composed of "Larkin domains" [Q of characteristic size much 



larger than the average inter-vortex spacing, and in a Bragg glass phase with algebraic decay 
of positional correlations. The fact that pinning should affect the lineshapes only weakly 
in such a state is expected on intuitive grounds because only the short-range order of the 
arrangement of the vortices is probed in /iSR experiments. However, these results disagree 
with experimental observations in the high-field limit, indicating that the system cannot have 
strong short-range crystalline order in this regime. In the absence of a detailed theory of the 
structure of the flux-line system in this regime, we propose that the positional correlations 
in this state resemble those in the liquid just above the solid-liquid "transition" . This state, 
however, differs from the high-temperature liquid in that the characteristic time scale for 
structural relaxation is much longer than that in the liquid. The correlations probed in the 
/iSR experiments are then naturally understood as being representative of the positional 
correlations in the equilibrated liquid. This proposal is consistent with the "pinned liquid" 
state postulated in the phase diagram suggested in Ref. [^]. This description should also 
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be appropriate if the high-field, low-temperature state of the system represents a new glassy 
phase that is separated from the high-temperature liquid by a line of true phase transitions. 
This is because the positional correlations in a structural glass are generally quite similar 
to those in the liquid near the glass transition. We present calculations of the second and 
third moments (the third moment is calculated using a specific decoupling approximation 
for correlations) of the field distribution in such a low-temperature disordered state and find 
that the results (summarized in Table I) are in good agreement with the experimental data 
on BSCCO. 

Several time scales are relevant to the analysis of transverse /iSR spectra in the mixed 
phase. A muon injected into the sample at time t = thermalizes due to interactions with 
the bulk material on timescales of order 10^^" s 0. This time scale is much shorter than the 
muon lifetime and typical precession times (~ 10^^—10^^ s). After thermalization, the muon 
diffuses in the sample. If the diffusivity of the muons is small, the typical distance traveled 
by a muon during its lifetime is much smaller than both the mean inter-vortex spacing 
and the penetration depth. In this limit, which is the one we consider, the polarization 
after time t is completely determined by the initial polarization and the distribution of 
local fields. (A complication that may be present in real experiments is that the muons 
may be preferentially absorbed at certain sites. We do not consider this possibility in the 
present work and assume that the muons are uniformly distributed.) Other time scales 
that enter the problem are the relaxation times of density fluctuations in the solid and 
liquid phases, which determine the time scales over which the field felt by a stationary 
muon varies. For the pure system, consistent with other work and with results obtained on 
YBa2Cu307__5 (YBCO) and other high-Tc superconductors, we assume that the relaxation 
time in the low-temperature (lattice) phase is very short relative to the inverse precession 
rate of the muons. (Experimentally, these time scales are found to be of order 10"^°s ^ 
close to the fiux-lattice melting transition.) In the fluid phase, the /xSR hnewidth would be 
expected to be vanishingly small due to motional narrowing if the vortex lines are highly 
mobile. However, if the vortex lines have very low mobility, such as might be expected in 



strongly pinned, entangled or glassy states, the second moment of n{B) can be related to 
the static structure factor of the liquid. We argue that the low mobility situation is realized 
at low temperatures in the high-field state of the pancake vortex system in the presence of 
random point pinning. At high temperatures and deep into the fluid regime, we argue that 
the vortex lines are highly mobile, so that the structural relaxation times of the vortex system 
are much shorter than the inverse precession rate of the muons. Even in this regime, there 
is a non- vanishing contribution to linewidths that comes from a space-averaged correlation 
function of pinning-induced density inhomogeneities in the liquid. We present a calculation 
of this contribution to the second moment of the field distribution function. Our results are 
representative of the hnewidths obtained experimentally deep into the liquid phase. 

The outline of this paper is as follows. In section |I|, we briefly discuss the calculation 



of the field distribution and its moments in a perfect Abrikosov lattice. In section |T|, 
we introduce our model of "pancake" vortices confined to the superconducting planes, and 
briefly describe the liquid-state and density functional methods used in our study. We also 
discuss the replica method we use to study the disordered liquid and its freezing transition 
in the presence of random point pinning. In section |V], we calculate the magnetic field 
distribution and its moments for a thermally broadened Abrikosov lattice. We use the 



m 



density functional description of the crystalline solid phase of the vortex system ||35| , |36 
our analysis. In section 0, we calculate these moments in a "frozen" liquid state of the 
vortex system. We provide calculations of the second and third moments to support our 
proposal that values of [Ai?^] consistent with Xab values of about I8OOA can be obtained if it 
is assumed that the low-temperature state of the vortex system at high fields is a glassy one 
characterized by short-ranged, liquid-like correlations and long relaxation times. The effects 
of weak disorder on the second moment of the field distribution in the high-temperature 
liquid state are discussed in section [V^, where we point out that non-vanishing linewidths 
can be obtained in this state by considering the contribution from the static inhomogeneities 



induced in the liquid by the quenched disorder. Section |V11| contains a summary of the main 



results obtained in this paper and a few concluding remarks. 
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II. /zSR SPECTRUM FOR AN IDEAL ABRIKOSOV LATTICE 

The time evolution of the transverse polarization (P+(r,t) = Px{r,t) + iPy{r,t)) of a 
muon in a locally inhomogenous magnetic field B(r), assumed to lie in the z direction (i.e. 



B(r) = zB{r)), is given by the Bloch-Torrey equation |44 

^^^^^ = -z7M5(r)P+(r, t) + DV'P^iv, t), (2.1) 

where D is the muon diffusivity, Px(r, t), Py{r, t) are the two components of the polarization 
transverse to the field and 7^ is the muon gyromagnetic ratio (= 8.516x10^ rad T^^ s^^)- 
Eq.( pyTD is to be solved subject to the initial condition P+(r, t = 0) = Pj^(r). If the initial 
distribution of muons is uniform then P+(r,t = 0) = P^ and if, in addition, the muon does 
not diffuse, the solution of Eq. (|2.1| ) is 



P+(r,t) = P°expH7^P(r)t]. (2.2) 

The space averaged distribution of the local magnetic induction is defined as 

n{B) = [5iB-B{r)], (2.3) 

where [^(r)] refers to the space average [l/V) Jy drA{r), V being the volume of the sample. 
It then follows that 

p(+)(t) = / n[ — \ exp{-iujt)diu, (2.4) 

7/. J'^ V7m/ 

where P^~^\t) is the spatial average of P^^\r,t) and 7^P = to. The measured P^^\t) 

can then be used to obtain, through an inverse Fourier transformation, the distribution of 

internal fields in the sample. 

The local magnetic field at a point r due to an array of straight vortex lines in the z 

direction is given by 

B{r)= fdr' b{\r-r' \)p{r'), (2.5) 
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where 6(r — r') denotes the field at point r due to a single vortex line at r' and p(r) is the fiux 
line density at point r (r is a two-dimensional vector perpendicular to the z-axis). Using the 
Fourier decomposition p(r') = (l/A) X]kPkexp(zk ■ r'), and 6(r) = (l/A) X]k^kexp(zk ■ r) 
(A is the cross sectional area of the sample) we get 

B{r) = {l/A)J2b^pkexp{{k-r). (2.6) 

k 

In an Abrikosov lattice, only terms involving reciprocal lattice vectors (k =G) survive in 
Eq.( p.6|) . For a perfect lattice with (5-function densities, pc = N for all G (A^ is the number 
of fiux lines) and the sum may be written as 

^W = |^E^Gexp(.G-r), (2.7) 

where A^$o = BqA and Bq is the magnetic induction within the sample, averaged over 
distances of order the penetration depth. $o is the quantum unit of fiux carried by each line 
and equals 2.07xl0~^^ T m^. In the London model, 

"o = ^^. (2.8) 

with 0^ = 0"^ + Gy, which leads to a magnetic induction that diverges logarithmically at 
the core of a vortex. This unphysical divergence can be eliminated by introducing a core 
"form factor" f{G) which represents the behavior near the core more accurately. This then 
gives 

Analytic expressions for f{G) can be obtained from solutions of the full Ginzburg-Landau 
equations at low inductions. A variational solution to these equations by Clem [^ yields 
the result f{G) = gKi(g), where g = ^/2^(G^ + A^^)^'^. Here, ^ is the coherence length 
of the superconductor and Ki{x) is a modified Bessel function. A less accurate but more 
commonly used representation of this core form factor replaces gKi{g) by the Gaussian form 
exp(— ,^^G'^/2). Recently, Yaouanc et al [^ have argued that it is essential to retain the full 
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form of the core form factor for a satisfactory treatment. Using results originally derived by 
Hao et al [^, they obtain 

«„.«„X:«^«^exp,G.r,. (2.10) 

where w'^{G) = 2{^^G^){l + b^)[l~2b{l-b)% b = Bo/B^2, and we have specialized their ex- 
pressions to the case in which there are no in-plane anisotropics. For the triangular Abrikosov 
lattice, the sum in Eq. (|2.1CI| ) is to be carried out over the reciprocal lattice generated by the 
vectors Gi = {An / ^/3ao)y and G2 = {2Tc/ao)[x + y/Vs]. Here the inter- vortex spacing oq is 
given by Oq = 2^q/\^Bq. For the calculations reported in this paper, we use the full form 
factor proposed by Yaouanc et al. in our description of lineshapes and linewidths in both 
solid and fluid phases. Also, since we work at inductions much less than Bc2, we set 6^ ~ 
and [1 - 26(1 - bf] ~ 1. 

The moments of the field distribution function n{B) can be obtained from Eq.( p.lO| ). 
The second moment of the distribution function [Ai?^] = J^^ n{B) {B — Bo)^dB is obtained 
as 

' '° °£ (1 + A^GT -''"S ^' ■ <"" 

where 

G^ = ^^[m^ + mn + n^], (2.12) 

OCLn 

m, n are integers and we have used A/ao ^ 1 in obtaining the final answer in Eq.( |2.11| ). The 
third moment, defined by [A5^] = J^^ n{B){B — Bq^cIB is obtained as 

rAR3l_R3 V- a{Gi)a{G2)cT{G3) 

^ ^" ° o„c&3.o (l + A^G?)(l + A2Gl)(l + A^Gi)' ^'-^'^ 

Gi+G2+G3=0 

where the cr{Gi), i = 1, 2, 3, represent the cutoff corrections and are given by 

a{G,)=w{Gi)K,{w{G,)). (2.14) 
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III. MODEL AND METHODS OF CALCULATION 

In extremely anisotropic, layered superconductors in their type-II regime in an external 
field applied perpendicular to the layer plane, the statistical mechanics of the interacting flux- 
line system reduces to the problem of interacting planar vortices. These "pancake" vortices 
are formed at the intersections of the flux lines with the planes on which the superconducting 
order parameter \l'(r) attains its maximum value. In the limit of infinite anisotropy (a 
reasonably accurate starting point for BSCCO in which the anisotropy factor 7 has been 
estimated to be higher than 100), the pancake vortices interact via a pair- wise potential that 
can be calculated within the framework of a London model ^^. This limit corresponds to 



assuming that the direct Josephson coupling between two adjacent superconducting layers 
is vanishingly small. However, vortices on different layers continue to be coupled through 
their electromagnetic interaction. Recent studies P,[l6| have shown that in BSCCO, the 
electromagnetic interaction between vortices dominates over the interaction generated by the 
Josephson couphng if the temperature is lower than about 0.8Tc(0). As the low-temperature 
behavior will be our primary concern in this paper, our assumption of infinite anisotropy 
can be expected to be quite reasonable for BSCCO. 

As we have argued earlier [^,Q, the statistical mechanics of the flux-line system in 
its fluid phase is most easily approached in this limit by considering the system of pancake 
vortices to be a classical, anisotropic liquid. The equilibrium density correlation functions of 
this classical liquid may be calculated by generalizing methods of liquid-state theory, and the 
freezing transition of this liquid may be studied using well known methods of classical density 
functional theory with the calculated liquid-state correlation functions as input. Such an 
approach makes it possible to characterize thermally broadened density distributions in the 
equilibrium solid at the freezing transition without using any ad hoc assumption. 

In this section, we first define the effective "Hamiltonian" of the vortex system in the limit 
of infinite anisotropy, and then provide brief introductions to the liquid-state and density 
functional methods used in our analysis. Finally, we discuss the techniques we use to study 
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the effects of disorder on liquid state correlations and the extension of density functional 
methods to deal with the problem of freezing in the presence of random pinning. 

A. Interaction between pancake vortices in the limit of infinite anisotropy 

If the coherence length of a layered superconductor in the direction perpendicular to 
the layers is of the order of the inter-layer spacing d, the appropriate generalization of the 
Ginzburg-Landau functional is a free-energy functional suggested by Lawrence and Doniach 
||37| . In the limit of infinite anisotropy (i.e. vanishing Josephson coupling between the 
superconducting layers), the Lawrence-Doniach functional in the "phase-only" (London) 
approximation reduces to a particularly simple form from which the interaction energy of a 
collection of pancake vortices can be readily derived [^ , ^ . The energy turns out to be a 
sum of pairwise interactions (this simplification occurs only in the limit of infinite anisotropy 
- the effective "Hamiltonian" of a system of pancake vortices can not be broken up into a 
sum of pairwise terms if the Josephson coupling is finite). In Fourier space, the inter- vortex 
interaction is given by 

.^. ^ ^ rA^[fci + (4/d^)sin^(M/2)] 

^^^ kl[l + \^kl + A{\y<P)sin\kJ/2)Y ^ ^ 

with r = pd^l/inX^ (A is the in-plane penetration depth - we have removed the subscript ab 

for ease of notation) and (3 = l/fc^T. In Eq.( p.lD , kj_ and k^ refer, as before, to directions in 

the plane of the layers and normal to them. Two vortices lying on the same layer repel each 

other with a potential that is a logarithmic function of their separation, whereas two vortices 

on different layers attract each other with a potential that is weaker than the intralayer 

potential by the factor d/X. The interlayer potential is logarithmic in the in-plane separation 

and falls off exponentially in the z-direction as e~"''^'^ where n is the layer separation. 

The Fourier transforms of the radial and 2;— axis components of the magnetic field due 

to a single pancake vortex at the origin are given by 

d<l>ok,k^ 



kl{l + X^kl + X^k^,)' 
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(3.2) 



and 

The calculation of local field distributions and their moments in this case proceeds in much 
the same way as outlined in the previous section. In this description, however, there is also 
a radial component associated with the total magnetic induction. (This component vanishes 
in the case of perfectly aligned straight vortex lines). In addition, the 2;-axis component of 
the magnetic induction is a weak function of z, even for a straight line of pancake vortices. 
If the inter-layer spacing is small, however, in comparison to the length scale over which the 
field varies (i.e. d ^ A), this variation can be neglected. 

B. Liquid state correlations and density functional theory for the pure vortex system 

A basic quantity that appears in a statistical description of the equilibrium proper- 
ties of a simple isotropic liquid is the pair distribution function g{\ r — r' |), defined by 
{J2iLi I^j^i '^('^■~i"i)'^(i^'~i^i)/P£)- Here, Fj denotes the location of the ith particle in the hquid 
and Pi is the average density. The pair correlation function h{r) defined by h{r) = g{r) — 1, 
and the static structure factor defined by S{k) = 1 +p£ J drhir) exp(ik ■ r) are closely related 
to g{r). The pair correlation function can be decomposed as 

h{r) = C{r) +pj dv'C{\ r - r' |)/i(r'), (3.4) 



which is the Ornstein-Zernike relation ^^ . The Fourier transform of the direct pair corre- 



lation function C{r) is related to the structure factor through S{k) = 1/(1 — piC{k)). 



These definitions can be readily generalized ||35| , |36[| to describe the equilibrium correla- 
tions in a layered liquid of pancake vortices. Since this system is highly anisotropic, the 
two-point correlation functions in position space have two arguments: the in-plane separa- 
tion p and the separation nd in the z-direction {n is an integer). The Fourier transforms 



of these functions depend separately on k^ and k^. As discussed in detail in Refs. PH] and 
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standard methods of liquid state theory, such as the Hypernetted Chain approxima- 
tion (HNC) [^ and the Rogers- Young closure scheme |5y], can be generalized to calculate 
these correlation functions of the layered vortex liquid. These correlation functions are then 
used as input to a density functional theory of freezing from which the freezing temperature 
and the time-averaged density distributions in the solid phase at the freezing transition are 
obtained. 

In the density functional theory, the free energy of an inhomogeneous configuration of the 
time-averaged density field is expressed in terms of fluid-phase correlation functions. Freez- 
ing into a crystalline structure occurs when the free energy of the state with the appropriate 
periodic density modulation equals the free energy of the uniform liquid. In the density 



functional theory of Ramakrishnan and Yussouff [Q, the grand-canonical free energy cost 
(with respect to the uniform liquid) of producing a time-averaged density inhomogeneity is 
expressed as a functional of the density p(r). The simplest such functional for an isotropic 
system is 






dr 



p(j) In ^ - 6p{r) 



- fdr fdr'C{\ r - r' \)Sp{r)Sp{r') + ... 



(3.5) 



where p(r) is the density at point r, Sp{r) = p(r) — pi, T is the temperature, pe is the liquid 
density and the ellipsis denotes higher order terms which are conventionally set to zero. 

In mean-field theory, density configurations which represent the equilibrium phase min- 
imize Eq. (|3.5| ), i.e. satisfy SAQ/6p{r) = 0. This condition gives 

-P(r) 



In 



JdrV{\r-r'\Mr')-p,]. 



(3.6) 



A periodic (crystalline) density configuration is a solution of the mean-field equations if the 
following self-consistency condition is satisfied: 



I + 1] + ^ Pg exp(zG ■ r) = exp 



CoV + Yl CgPg exp(zG ■ r) 

G^O 



(3.7) 



where the Gs are the reciprocal lattice vectors of the structure to which the liquid freezes, 
Pgs are the Fourier components of the density field with wave vector G and we have defined 
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C = piC. The quantity t] is the fractional volume change on freezing from the liquid. Note 
that the uniform liquid, for which pc = for all G 7^ and 77 = 0, is always a local minimum 
of the free energy. However, because of the non-linearity of the self-consistency condition, 
periodic density- wave solutions which could be absolute minima of the free energy of Eq. ( p. 51) 
may appear as the correlations in the liquid increase. The pcs which minimize the density 
functional represent the (mean-field) density distribution in the crystalline solid. Since the 
density distribution near a lattice point is not restricted to be Gaussian in a fully self- 
consistent theory, the thermal broadening of the peaks of the density distribution includes 
at least some of the anharmonic effects that are expected to be present near the melting 
transition. Our treatment, therefore, goes beyond the harmonic approximation made in 



several existing studies [|7|,|8|, |29|j30|| of /iSR spectrum in the mixed phase. 

The density functional formalism outlined above can be generalized in a straightforward 
manner [^,0] to a layered system of pancake vortices. The reciprocal lattice vectors that 
characterize the Abrikosov lattice are (G^,^^ = 0) where the G^s are two dimensional 
reciprocal lattice vectors defining a triangular lattice. 

C. Replica liquid state theory and density functional theory in the presence of 

pinning disorder 



In earlier work |^ , two of us developed a replica theory of the correlations of a vortex 
liquid in the presence of random point pinning and derived a replicated free-energy functional 
analogous to the Ramakrishnan-Yussouff functional to describe the freezing of this disordered 
liquid. We briefly summarize our methods here. Our analysis is based on the replica method 
||5T| applied to a system of classical particles interacting via the Hamiltonian 

H = Hk^neuc + ^ E ^(1 r.. - r, I) + E Mrt), (3.8) 

where V{r) is a two-body interaction potential between the particles and Vd{r) is a quenched, 
random, one-body potential. We assume that Vrf(ri) is drawn from a Gaussian distribution 



-| „ -I n n i\ i\ 

'^- uK,J^<-M-^rTT.T.T.T.y'''i\-t--U))- (3-9) 

"^■J J ^^bJ- „=i 3^1 i=i ,=i 



of zero mean and short ranged correlations: [Vrf(r)Vd(r')] = K{\ r — r' |), with [■ ■ •] denoting 
an average over the disorder. Using [In Z] = hni„_^o[(^"' ~ !)/'"■]) o^is obtains, prior to taking 
the n —>■ hmit, a rephcated and disorder averaged partition function of the form 

T n n N N 

Here a,/3 are rephca indices and V^^d r" — r^ |) = V{\ r" — r^ \)6ai3 — PK{\ r" — r^ |). 

Our approach to the problem defined by Eq. ( p. 91) begins by recognizing that it resembles 
the partition function of a system of n "species" of particles, each labeled by an appropriate 
replica index. These particles interact via a two-body interaction which depends both on 
particle coordinates (rj,rj) and replica indices {a, (3). This system of n species of particles 
can be treated in liquid state theory by considering it to be a ra-component mixture and 
taking the n — > limit in the Ornstein-Zernike equations governing the properties of the 
mixture. In the n -^ limit, assuming replica symmetry, we get the equations 

^'~ [l-peCWik)+PeC(^KkW ' ^ ^ 

and 

These equations are written in terms of the pair correlation functions h°'^ and the direct 
correlation functions C"'^ of the replicated system. The assumption of replica symmetry 
implies that C"'^ = C^^)6o,p + C^^\l - 6o,p) and /i"^ = h^^^6o,p + h^^\l - 6o,p). The physical 
interpretation of these correlation functions is as follows. The function /i^^^ describes the 
disorder-averaged equal-time (equilibrium) correlation of fluctuations of the local density, 
and h^'^^ represents the disorder-averaged correlation of disorder-induced deviations of the 
time-averaged local density from its average value p£. 

In Ref. [^, we calculated the functions h'^^\p,nd) and h^'^\p,nd) for the layered vortex 



system using the HNC closure approximation. To estimate the inter-replica interaction 
PK{p,nd), we assumed that the principal source of disorder is atomic scale pinning centers 



such as oxygen defects ||38| which act to reduce Tc locally. This yields 
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PK{p,nd) ~ r'exp(-pVe')5n,o. (3.12) 

Here ^ ~ ISA is the coherence length in the a6-plane, and F' ^ 10~^r^ for point pinning 
of strength dr^H^/STi, where vq is an atomic scale length parameter and we have assumed 
defect densities of the order of 10^°/cm'^. In Figs. ^ and ^, we show typical examples of 
the calculated correlation functions g^^^ and g^"^^ obtained using this theory at two different 
values of T at fixed induction Bq. 

The calculated correlation functions C^^\r) and C^'^\r) can then be used as input into 
an appropriately generalized version [|^ of the density functional theory for a mixture of n 
species of particles. This leads to the following density functional in the n —>■ limit: 



ar 



ksT 



p(r) In ^ - 6p{r) 



-ljdrJdr'[C^'\\r-v'\)-C^'\r-r'\)Mv)-peMr')-pe].... (3.13) 

Here we have assumed that the density field is the same in all the replicas (p"(r) = p(r) for 
all a). 

Generalization of this free energy to a layered system of pancake vortices is straightfor- 
ward. In Ref . ^l| , we assumed that if the pinning disorder is weak, then the "solid" phase of 
the vortex system is described by a global minimum of this free energy characterized by non- 
zero values of the crystalline order parameters pc,. Since we are interested in a mean- field 
description of a first-order transition (the fiux-lattice melting transition is experimentally 
found to be first order if the pinning disorder is weak), this assumption is justified as long 
as the low-temperature state exhibits strong local crystalline order (specifically, as long as 
the translational correlation length in the low-temperature phase is much longer than the 
range of positional correlations in the liquid just above the freezing transition). Results of a 
full calculation of the phase-boundary obtained through this method were presented in Ref. 
| 4T| . Our principal result was that the liquid-solid phase boundary in the {B — T) plane in 



the absence of disorder is only mildly affected by the presence of weak disorder, the sup- 
pression of this phase boundary to lower temperatures being larger at higher fields. These 
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predictions have been confirmed in recent experimental studies [p3] of this phase boundary 



in electron-irradiated samples of BSCCO. In the present work, we use the results obtained in 



Ref. [^ for the density distribution in the "nearly crystalline" phase to analyze the effects 
of weak pinning disorder on /xSR lineshapes in the low-temperature, low-field state of the 
vortex system. 

IV. /iSR SPECTRUM IN THE CRYSTALLINE PHASE OF BSCCO 

In this section, we first consider the crystalline (Abrikosov lattice) phase of the vortex 
system in the absence of pinning disorder. If the typical time scales for phonon-like density 
fiuctuations in this state are much shorter than the characteristic time scale of precession of 
the muon spin (~ {j^B)^^), the muons see a broadened (equilibrium) density distribution 
of fiux lines. The time-averaged local magnetic field B{r) is then obtained through the 
convolution relation 

Bir) = Jdv'b{\r-v'\)piv'), (4.1) 

which differs from Eq. (|2.5|) in that p(r) now represents the time- averaged density distribution. 
The Fourier components of the periodic time-averaged density field are given by Npo, where 

rV 



j-V 

Pg = drp{r)exp{-iG -r). (4.2) 



Here V denotes the area of the unit cell. In thermal equilibrium, this time average is equal 
to a thermal average. 



We calculate pa by solving, in a one-order-parameter approximation |^, the self- 
consistent equations for the Fourier components of the periodic density field at freezing. 
The input direct correlation function, C{k±, kz), is taken from a solution of the liquid-state 



equations with a Rogers- Young closure scheme [^,Q. We use parameters appropriate to 
BSCCO i.e. A(T = 0) = 1500 — I8OOA and d = 15A and assume a two-fiuid temperature 
dependence of A with Tc(0) = 85 K. As mentioned earlier, the relevant input correlation 
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function for the freezing of the vortex hquid into an Abrikosov lattice is C{G±,Gz = 0). 
In Fig. 1^, we plot the Fourier components pc of the periodic density field (triangles) at the 
freezing temperature as functions of G. Note the fast decay of pa at large G and the nearly 
Gaussian envelope of the decay. In Fig. |^, we display the magnetic field distribution for the 
thermally broadened Abrikosov lattice at freezing for values of the magnetic induction of 
(a) Bq = 0.3 T and (b) Bq = 0.05T, incorporating corrections due to the finite core size (see 
Eq.( p.lOD ). The field distributions obtained in this theory are notably narrower and more 



symmetric than the corresponding distributions calculated for an almost perfect Abrikosov 
lattice shown in the same figure for comparison. This feature is seen in experiments on 

Bscco igm. 

Although the density functional description is in principle applicable even deep within 
the solid phase, it is unclear how to parametrize the appropriate input correlation functions 
in this case. While it is well-known that the existence of the liquid-solid transition cannot 
be inferred from usual liquid state methods such as the ones used by us (the correlations 
show no singular behaviour), using these correlation functions in a density functional de- 
scription of the solid much below its crystallization temperature is questionable. The issue 
of how to describe the solid at temperatures far from freezing is relevant in determining the 
temperature dependence of the crystalline order parameters pa, which in turn determines 
the temperature dependence of the second moment of the field distribution function. 

In this calculation, we extract the temperature dependence of the pc's by assuming that 
the square of the width of the equilibrium density distribution at a lattice site is proportional 
to the temperature T at all temperatures lower than the melting temperature. The constant 
of proportionality is obtained by demanding that the solid at its melting temperature has 
the same Lindemann parameter as that obtained from our self-consistent density functional 
calculation at freezing. A linear temperature dependence of the mean square displacement 
also occurs in a harmonic treatment (with temperature-independent elastic constants) of the 
fluctuations of a vortex line about its equilibrium position at a lattice site. However, our 
treatment differs from a harmonic one in that the value of the mean square displacement 

22 



at the melting temperature is obtained from our density functional calculation which, as 
argued above, is non-perturbative and includes possible anharmonic effects. 

The thermal broadening of the density distribution in the crystalline phase is cal- 
culated in the following way. Our density functional calculations p5|,|36[] show that at 



fields larger than about 50 mT, the Lindemann parameter J {u^) /a^ at the freezing tem- 
perature is nearly constant at a value ~ 0.2. Assuming {u'^{T)) /a^ = xT, and using 
the value of the melting temperature (~ 18K) obtained from our calculation, we ob- 
tain X — ^/4:50K~^. If we assume that the density distribution in the crystalline state 
is represented by a sum of Gaussian density profiles centered at the lattice sites, i.e. 
p(r) = Z]j /(i" ~ R-i)? with /(r) = a/7rexp(— ar^), the assumed form of (m^) implies that 
a(T) = l/{alxT)- The time-averaged Fourier components of the density field are then given 
by Pg = exp(— G'^aQxT/4) = exp(— 47r^(?7i^ -|- n^ + mn)xT/3). The second moment of the 
local field distribution is calculated using 

The temperature dependence of the second moment comes from two sources - the (assumed 
two-fluid) temperature dependence of A on T, as well as the temperature dependence of 
the pgs. Fig. |^ shows the temperature dependence of the second moment as calculated in 
this approximation. Note the curvature of the plot of (AB)"^ vs. T. This trend is opposite 
to that obtained if one assumes a two-fluid form for the T-dependence of A and ignores 
the corrections arising from the temperature-dependence of the broadening of the density 



distribution. A similar trend is seen in some experiments |jT3| at not too high fields {Bq < 1.5 
T). While the use of Debye- Waller-like broadening factors to rationalize the non- two- fluid 
temperature dependence of [Ai?^] is not new (see |T3|J30[]), in our calculation, these factors 



are not obtained as fitting parameters but follow from the general structure of our theory. 

In Fig. 1^, we show the temperature dependence of the third moment, as obtained in this 
calculation using 
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[Ar3i = r3 V '^i.Gi)a{G2)(j{G^)pG^pG^pG^ 

Gi+G2+G3=0 

where <7{Gi), i = 1,2,3 represent the cutoff correction described in section |T| (see Eq.( |2.f4| )). 



The resuhs of our rephca treatment |^ of the effects of weak pinning disorder on the 
correlations and freezing of the vortex system indicate that the results obtained above for 
the pure system would remain essentially unchanged if the low-temperature "solid" state in 



the presence of pinning exhibits strong local crystalline order. As discussed in section |lTl| C 
above, the freezing of the liquid into such a state may be described by a density functional 
(see Eq.( p.l3D ) in which the quantity [C'^^^r) — C^'^\r)] plays the same role as that of the 



direct correlation function C(r) in the density functional theory of freezing of the pure 
system. Our replica liquid state theory |^ shows that if the pinning disorder is weak, 



then the function G^^\r) is nearly identical to the function C(r) for the pure system, and 
G^'^\r) is much smaller, so that C^^\r) — C*^^^(r) ^ G{r). This, in turn, implies that the pcs 
that describe the (local) density distribution in the low-temperature state in the presence of 
pinning are nearly the same as those obtained for the pure system. Since the /xSR spectrum 
in the low-temperature state is determined by the values of the pcs (see Eqs. ( |4.3| ) and 
([4.4D), it would be nearly unaffected by pinning disorder as long as the approximation of 
describing the local translational order in the low-temperature state by the pcs remains 
valid. This is expected because, as discussed earlier, pSR spectroscopy probes the local 
translational order of the vortex system. 

Before we conclude this section, we comment briefly on the liquid-solid phase boundary 
obtained in our calculation in relation to the experimental data. Our calculations, using a 
Rogers- Young closure scheme and a one-parameter density functional theory, indicate that 
the liquid-solid phase boundary in the pure system at high fields is practically independent 
of field and occurs at T^j ~ 18K. This value is slightly lower than the values obtained in 
experiments @,0]. In our calculations, we used values of Xab = 1500— I8OOA. Assuming larger 
values of Xab would act to shift the phase boundary to still lower temperatures. However, 



as we have argued earlier |]36|, the principal source of error in our estimation of the phase 
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boundary arises from the fact that our hquid state approximations underestimate correlations 
in the hquid. A more accurate treatment of hquid state correlations is expected to shift 
the phase boundary to higher temperatures. Also, a comparison of our results for the pure 
system at high fields with the experimental data may not be particularly meaningful because, 
as discussed in section | above, the freezing observed in experiments at high fields is probably 
to a glassy state rather than to a crystalline one. 

V. /iSR SPECTRUM IN A "FROZEN" LIQUID 

We now discuss the calculation of the moments of n{B) in the liquid phase of the vortex 
system. If the time scale of structural relaxation in the liquid (i.e. the time scale over 
which the local density becomes homogeneous in the pure system, or nearly homogeneous 
if weak pinning is present) is much smaller than the inverse muon precession rate, the 
density distribution seen by the muons should be nearly homogeneous and the line widths 
vanishingly small. In the opposite limit, in which structural relaxation in the liquid occurs 
over times that are much longer than the inverse muon precession rate and the muon lifetime, 
some quantitative progress may be made in calculating /xSR lineshapes and linewidths. We 
propose that such a limit is appropriate for a description of the low-temperature, high-field 
state of the fiux line ensemble in BSCCO. Such glassy behavior may have several origins, as 
discussed in the Introduction. Irrespective of the precise nature of this state, it is a reasonable 
guess that it exhibits short range positional correlations but no long range crystalline order. 
A natural assumption then, is that positional correlations in this state should resemble 
those in a liquid. As mentioned in the Introduction, this proposal is consistent with existing 
neutron scattering [^] and other experimental data. It is also physically sensible, as it 



provides a way by which some degree of local correlations can be built into the glassy state, 
while ensuring that crystalline long-range order is destroyed, as required by the Larkin- 
Ovchinnikov argument |^2[. In this section, we explore the consequences of this proposal. 
In this (adiabatic) limit, the muons essentially see a "snapshot" of the liquid (i.e. a 
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typical liquid-like configuration of the pancake vortices). The second moment of the field 
distribution function is then given by the space- averaged two-point correlation function of the 
density in the configuration seen by the muons. If the system is self-averaging (we assume 
that it is so), then this correlation function is nothing but the time-averaged (equilibrium) 
two-point density correlation function of the liquid. Therefore, the linewidth in this limit 
can be calculated in terms of the static structure factor S{k) of the liquid, which is defined 



as 



5(k) = ^(PkP-k) = 1 + ^(Eexp(zk ■ r,,)), (5.1) 

where Tjj = Fj — Tj (r^ is the position vector of the ith pancake vortex), A^ is the total 
number of vortices, the sum runs over all the vortices, and the brackets (■ ■ ■) denote a 
thermal average. The second moment of the field distribution function is given by 

[AB^] = [ABl] + [ABl]. (5.2) 

The brackets [■ • ■] denote a space average for a typical liquid-like configuration of the vortex 
system. Precisely which moment ([Ai?^] or [Ai?^^]) is measured depends on the geometry 
in which the experiments are done. Our results are illustrated here for the transverse /xSR 
geometry, in which the muons are injected with their polarizations in the plane of the 
superconducting layers and thus perpendicular to the external field. In this geometry, the 
measured quantity is [Ai?^]. 

The equivalence of space and thermal averages then leads to the following expression for 
the second moment of the distribution of Bz'. 

[ABl] = ^T^ J d'k^dk,bl{k^,k,)S{k^,h), (5.3) 

where pi = Bq/^q denotes the (areal) density of vortices, and the quantity b^ is defined 
in Eq.( |3.3| ). A similar expression for [Ai?^] may be obtained by replacing b^ by b±_ (see 
Eq.(|3.2|)) in the above equation. 

Before analyzing the relevance of short-range correlations in the computation of moments 
of n{B) in the frozen liquid, we discuss two simple limiting cases. First consider the limit in 
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which the hquid has no correlations, so that S{k±, kz) = 1, and the second moment of the 
distribution of the ^-component of the field is given by 



,2, _ Pi H<i' r ,. ,,. r * 



For B = 0.5T, d = 15A and A = 1500A, we obtain [AS^] = 1830 G^ whereas assuming 
A = I8OOA gives [ABl] = 1059 G^. 

Now consider the case of a completely uncorrelated arrangement of rigid flux lines. The 
second moment of the field distribution is then given by 

[ABS = 7^/7^4Lt.A-, = ^. (5.5) 



(27r)2y {1 + X^kiy "- 47rA2 

Assuming, as before, Bq = 0.5T and A = 1500A, we obtain [AS^] = 3.66x10^ G^ (2.54 
xlO^ G^ for A = 1800 A). These numbers are to be compared with the value of [Ai?^] (= 
3140 G^ for A = 1500A, 1514 G^ for A = I8OOA) for a perfect Abrikosov lattice. 

These simple results, first derived by Brandt |5^, bring out an important point: the 
destruction of correlations in the a6-plane (while keeping the vortices perfectly correlated 
in the c-direction) greatly increases the linewidth, whereas the destruction of correlations 
in the c-direction causes the linewidth to become smaller. This observation suggests that 
a likely candidate for a structure that exhibits anomalously narrow /iSR linewidths would 
be one with at least local (see below) positional correlations in the layer planes and little 
correlation across the layers. Our proposal for the structure of the "glassy" state of the 
vortex system is consistent with this expectation: our liquid state calculations ||35|j36[] show 



that the equilibrium liquid just above freezing exhibits pronounced short-range correlations 
in the layer plane, but the correlation length in the c-direction remains small (about 15 - 
20 layer spacings). 

We now describe a simple calculation that illustrates that the presence of short-range 
a6-plane correlations in a frozen liquid of rigid flux lines is sufficient to bring down the 
the linewidth to values comparable to that obtained for a perfect Abrikosov lattice. The 
inclusion of local correlations in the calculation of the moments of the field distribution 
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requires the specification of the hquid structure factor S{kj_). At the simplest level, S{kj_) 
can be approximated as 

S{k^) = {0<k^<km-e), 

= Dm [km — e < k± < km + CJ, 

= 1 {km + e < k± < oo). (5.6) 

This approximation retains all the basic qualitative features expected of a fiuid structure 
factor - it rises from a value of zero near the origin of A;— space to a peak at a fc— space 
distance of the order of 27r times the inverse of the mean interparticle spacing, and becomes 
unity at sufficiently large k. With this approximation, it is easy to obtain 

Assuming 2eSm/km ~ 0.5, so that [AB^] ~ 5<l>o/(27rA^A;^), A;^ ^ An^B/^o, and for suffi- 
ciently large fields (so that (1 + A^A;^) ~ A^A;^), we obtain 

[Ai?.1=^^^4xlO-^||, (5.8) 

which is quite close to the result, [A_B^] ~ 3.71 x 10^^^, for the linewidth in a perfect 
Abrikosov lattice. This admittedly crude calculation thus brings out the important point 
that the /iSR linewidth is sensitive only to the local order. 

We now proceed to describe our detailed calculation of the linewidths in a frozen liquid- 
like state. In terms of the direct correlation function, Eq.( |5.3| ) can be written as 

f^^'] = f 12^0^/0 ^^^^^ loo "^^^ (1 + X'k{ + X^klY (1 - p,C{k^, A:,)) ' ^^-^^ 

where C{k_j_,kz) is the liquid direct correlation function and we have used the relation 
'S'(A;j_, kz) = 1/(1 — PeC{kj_, k^)). The direct correlation function C{kj_, k^) is obtained from 
our earlier calculation I^S],^ of the equal-time correlations in a layered vortex liquid. To 
represent the core form factors (in order to reduce the complexity of the expressions, we do 
not explicitly show these form factors in the equations derived below), we model the core of 
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a pancake vortex as a cylinder of radial dimension ^ab — ISA and height ^c — 3A. We use 
the decomposition 

C{p,nd) = Cs{p,nd) + CL{p,nd), (5.10) 

where the "long-range" part of C is given by CL{p,nd) = —f3V{p,nd) {V{p,nd), the inter- 
vortex interaction in position space, is the Fourier transform of V^(k) defined in Eq.( p.l| )). 
Asymptotically, C{p,nd) tends to the value —(3V{p,nd), which we have called Cl(p, nrf), 
so that Cs{p,nd) tends to zero. This "short-range" part is assumed to have the form 
Cs{p,nd) = Cs{p,n = 0)6n,o (see [ ^ , ^ ) . The Fourier transform, Cs{k±), of the func- 
tion Cs{p,n = 0) is obtained by numerically solving the self-consistent equations of liquid 
state theory. Approximating sin^(fc^(i/2) appearing in Eq. (|3.1|) by k'^d'^/A, the integral over 
kz in Eq. (|5.9| ) can be performed to yield 






AtxTX^ Jo ^ ^yi X \ i + x^k 



' ' (5.11) 



,2 ' 



where 



X = kl{l + \^kl){l - p,Cs{k^)) + p.TX^kl, 

Y = \^kl{l-p,C,{k^))+p,V\\ (5.12) 

These equations and the numerically obtained Cs{k_[_) are used to evaluate the second mo- 
ment [A52]. 

The third moment of the field distribution function in the liquid phase can be obtained 
via the liquid state methods we have discussed here, together with a simple decoupling 
approximation. Using 

\^Bl] = ^^ E (P(ki)p(k2)p(-ki-k2))6.(ki)6,(k2)6,(-ki-k2), (5.13) 

^ ki,k2,k3 

what is required to calculate this moment is the triplet structure factor 5'*^^)(ki, k2) defined 
by 

S(=^)(ki,k2) = l(p(ki)p(k2)p(-ki -k2)). (5.14) 
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In terms of the triplet direct correlation function C^(ki, k2), this assumes the form 

S(ki)S(k2)S(-ki - ka) 

The triplet correlation function is typically small in magnitude in regular three dimensional 
solids. If it is assumed to be zero, (the so-called convolution relation [^), the integral can 
be written completely in terms of the the static structure factor. We thus obtain 



[ABl] = ^-i- J dk,^ J dh, J dk2± J dk2. 



d {2ny 
S{k^)S{k2)S{-k, - k2)6.(ki)fe,(k2)6,(-ki - ka). (5.16) 

This expression is defined in terms of the two-particle structure factor S{k±, k^) which, 
as discussed above, may be obtained from the function Cs{k±) calculated in our liquid 
state theory. Exploiting the structure of this expression, and after considerable algebra, the 
6-dimensional integral is reduced to a 3-dimensional one which is evaluated numerically. 

Our results for [Ai?^] and [ASf] for values of Xab = 1500A and I8OOA are summarized 
in Table I. These results are obtained using the Cs{k±) calculated just above the melting 
transition. The core form factors, omitted for simplicity in the expressions derived above, 
are present in the full calculation. While our liquid-state calculations are performed for 
A = I500A, the calculated Cs{k±) can be used for the case Xab = I8OOA in the high field limit 
provided the melting temperature is scaled as Tm oc 1/A^^. This is a simple consequence 
of the crossover to quasi-two-dimensional behavior at high fields, in which the dominant 
interaction is the (scale-invariant) logarithmic coupling between vortices on the same plane, 
which is perturbed only weakly by the presence of vortices on other planes. Comparison 
with experimental data Jl^ shows that our proposal agrees quite well with the results of 



experiments. 

Results for [AS^] can be obtained from 



i + \^e. 



(5.17) 



' A2 

Here X and Y are as defined in Eq. (|5.12|) . The value of this moment is comparable to [Ai?^] 
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Our ansatz that the vortices remain stationary during the hfetime (or precession time if 
it is smaller) of the muons is expected to be strictly valid only at very low temperatures - 
at higher temperatures we would expect the effects of thermal broadening of the (now inho- 
mogeneous) density distributions to reduce linewidths further in analogy with our results for 
the crystalline phase in the pure system. There seems to be no simple way of incorporating 
such thermal effects into the theory described here (see, however, the discussion in section 
Vn] below). A theoretical investigation of the temperature dependence of the linewidth in 



this disordered phase would be interesting. It would also be worthwhile to carry out a more 
accurate calculation of the third moment of the field distribution by going beyond the simple 
decoupling approximation used here. 

In the calculations described above, we used the equilibrium correlation function of the 
pure liquid at a temperature just above the crystallization temperature Tm to determine the 
moments of the field distribution. Since the glassy phase we are interested in is expected to 
arise as a consequence of the presence of pinning disorder, one should, in principle, use the 
correlation function of the disordered liquid near its freezing (glass transition) temperature 
in these calculations. Due to reasons mentioned below, we expect that the results of such a 
calculation would be essentially the same as those described above. Our replica calculation 
pl| of liquid state correlations in the presence of pinning disorder shows that the equal-time 
(equilibrium) density correlations in the liquid are basically unaffected by weak disorder. 
Therefore, our approximation of replacing the static structure factor of the disordered liquid 
by that of the pure one is expected to be quite accurate. We also find that the value of the 
linewidth depends relatively weakly on the temperature at which the liquid-state correlations 
are calculated - a change in the temperature from Tm to STa/ leads to a decrease of about 
20% in the linewidth [Ai?^]. This observation indicates that a small difference between Tm 
and the actual glass transition temperature of the disordered liquid would not cause any 
significant change in the final results. 
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VI. fiSn LINEWIDTH IN THE DISORDERED LIQUID 

The time scales for structural relaxation in the high-temperature liquid phase are ex- 
pected to be short. The second moment of the field distribution function in this phase of 
the pure system would be vanishingly small if these time scales are fast compared to the 
inverse precession rate of the muons. Even in this limit, however, measurable line-widths can 
obtain in the presence of quenched pinning disorder. The source of this broadening is the 
off-diagonal (in replica space) density correlation function. Off-diagonal correlations, which 
are absent in the pure system, become non-zero in the presence of disorder. These correla- 
tions arise due to the disorder-induced inhomogeneity of the time-averaged local density. In 
a replica symmetric description, the off-diagonal radial distribution function is defined as 

g('\r) ^ ((p"(0)p^^"(r)))/p| = [(p(0))(p(r))]/p^ (6.1) 

Here the ps represent local number densities, a and P are replica indices, ((■ ■ ■)) represents an 
average over the canonical distribution function of the replicated system and (■■■), a thermal 
average for the original disordered system prior to averaging over disorder. The square 
brackets [■ ■ ■] denote an average over the probability distribution governing the disorder. 

If the relaxation times of the liquid are short, then the muons see the time-averaged 
(equilibrium) desnsity distribution which is inhomogeneous if pinning disorder is present. 
In this case, the second moment of the field distribution function can be obtained from 
the correlation function /i*-^^ = g^"^' — 1 that describes the correlation of the deviations of 
the time-averaged local density from its space-averaged value pe. Specifically, the second 
moment of the distribution of the 2;-component of the local magnetic field is given by 

[^Bl] = ^-^Jdk^JdKblS'^'\k^,h), (6.2) 

where S'*^^^(k), the off-diagonal structure factor defined as [{Spk){5p_k)]/N, is given by (see 
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Using this and the expression, Eq.(|3.3|), for bz, Eq.(|6.2D can be written as 

r . p2i _ P^ 1 f ,, [ ,, ^0^^ PeC^^Kk^,h) . . 

Separating the diagonal direct correlation function C^^-* as before into a short-range and a 
long-range part and evaluating the integral over kz, we get 

where P and Q are given by 

P = (1 + \'kl){l - p,CP{k^) + p,C^'\k^)) + p,TX', 

Q = A2(1 - piCi'\k^) + p,C^'\k^)) + ^. (6.6) 

In deriving the above equation, we have used the approximation C^'^\p,nd) = C^'^\p)5nfi, 
so that its Fourier transform C^'^\k±,kz) does not depend on kz- This approximation is 
justified [^ because the interaction between different replicas (Eq.( p.l2"D ) does not couple 



vortices on different layers. 

Evaluating the integral in Eq.(|6.5|) numerically, using the values for C'-^-* and C^^-* ob- 



tained in Ref. |^ through a self-consistent solution of the replicated liquid state equations, 
we obtain values for [AS^] in the range 1 — lOG^. At Bq = 0.3 T, just above the melting 
transition, our calculations yield [Ai?^] = 4.76*^. The smallness of this number reflects 
the numerically small value of C^'^\k±). Further, this quantity vanishes as the off-diagonal 
coupling (equivalently, the disorder strength) is reduced to zero, leading to the expected 
vanishing of the linewidth in the high-temperature liquid phase of the pure system. The 
relatively strong dependence of C^^-* on the disorder strength suggests that the linewidths ob- 
tained in the liquid phase at high temperatures should increase as the sample is progressively 
made more disordered. 

VII. SUMMARY AND DISCUSSION 

In this paper we have presented a theory for /iSR spectra in anisotropic, layered super- 
conductors which incorporates thermal effects in both the solid and the liquid phase, and the 
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effects of weak pinning disorder. In particular, we have computed the effects of thermal fluc- 
tuations and disorder in determining the /iSR linewidths and lineshapes in the mixed phase 
of the extensively studied high-temperature superconductor BSCCO. In the crystalline solid 
phase, our density functional description of the time-averaged density distribution offers a 
nonperturbative way of assessing the effects of thermal broadening. Our theory of /iSR spec- 
tra in the liquid phase represents, to the best of our knowledge, the first attempt to calculate 
linewidths in this phase beyond the simplest assumption of uncorrelated lines or vortices. 
Our analysis suggests that the experimentally obtained linewidths in the low-temperature, 
high-field state in BSCCO may be understood if one assumes that vortices in this regime 
are in a frozen state that resembles a liquid in its local correlations. We have also quantified 
the effects of pinning disorder on the /iSR linewidth in the high-temperature vortex liquid. 
The results of our calculations are in qualitative (quantititative in some cases) agreement 
with those obtained in experiments on BSCCO. 

An implicit assumption we have made is that the flux-lattice can be treated classically, 
i.e. the characteristic temperature below which quantum effects are non-negligible lies far 
below the equilibrium melting temperature. This assumption has been questioned recently 
by Bulaevskii it et al. |^B| who argue that quantum effects are non-negligible except at 
temperatures close to the superconducting transition temperature Tc(0). Their argument 
rests on the assumption that the characteristic energy scale for quantum effects is set by 
the superconducting gap and therefore is ~ ksTci^O). While much is still to be understood 
about the precise role of quantum effects in the mixed phase, we believe that the success 
of the density functional approach in describing flux-lattice melting (which assumes that 
a classical liquid state description can be used in treating the liquid) serves as a counter 
example to the assertion that quantum effects must be included in any theory for the mixed 
phase. As we have emphasized, our theory yields a value of the freezing temperature that is 
lower than the experimentally observed value - quantum effects would surely act to depress 
this value further. 

Our results for the effects of pinning disorder on the /iSR spectrum are based on the 
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replica liquid state theory developed in Ref. |41] where only replica symmetric solutions of 



the liquid state equations were considered. A recent study [Q of the same set of equations 
for simple isotropic liquids shows that these equations exhibit a replica symmetry break- 
ing transition which may be interpreted as a glass transition. This work suggests that the 
low-temperature, high-field state of BSCCO may correspond to a replica symmetry broken 
solution of the liquid state equations derived in Ref. |^ . In Ref. [|^ , the two-point corre- 



lation function of the "frozen" local density in the glassy phase is found to be quite similar 
to the equal-time (equilibrium) correlation function of the local density in the liquid phase. 
This observation supports the ansatz we have used in this paper for the structure of the 
low-temperature, high-field state of BSCCO. 
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TABLES 
TABLE L Summary of our results for [Ai?^] and [Ai?^] ^^ a slowly relaxing liquid of pancake 
vortices just above its equilibrium freezing temperature Tm- The results shown are for Bq = 0.3 
T, and two values (1500 A and 1800 A) of the a6-plane penetration depth A. According to our 
ansatz (see text) for the structure of the low-temperature disordered phase found in the high-field 
mixed-phase regime of BSCCO, these calculated values should correspond to the experimentally 
observed data in this phase at low temperatures. Experimental /xSR spectra for BSCCO crystals 
at 5 K show, for Bq = 0.3, 0.4 and 1.5T respectively, squared linewidth values of [Ai?^] = 90 it 4, 
80 ± 5, and 52 ± 8 G^ and third moments of [AS^j ^ ;lOOO ± 80, 500 ± 100 and ± 250 G^ (see 



Ref. |13|,|0| 



A(l) 


[Ai?2] (G2) 


[ABl] (G3) 


T]\j (calculated) 


1500 


175.3 


4022.0 


~ 18 K 


1800 


111.9 


1117.5 


~ 12 K 
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FIG. 1. Schematic H — T phase diagram of a type-II superconductor in mean- field theory, 
ihustrating Meissner, mixed and normal phases separated by continuous phase transitions. 
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FIG. 2. Schematic H — T phase diagram for pure BSCCO illustrating the effects of thermal 
fluctuations on the mixed phase. Note the large region of "flux liquid" phase separated from a 
low-temperature crystalline phase by a first-order liquid-solid phase boundary. A very narrow 
region of liquid phase expected just above the the Hd line is not shown. The mean-field Hc2 
curve represents a crossover between a regime with strong amplitude fluctuations and one in which 
amplitude fluctuations are very weak. 
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FIG. 3. Schematic phase diagram for BSCCO incorporating effects of thermal fluctuations and 
quenched disorder due to random point pinning. The phases shown are the "Bragg Glass" phase (I), 
the "pinned liquid" or "vortex glass" phase (II), the disordered liquid phase (III) and the normal 
phase (IV). The existence of the Bragg glass and vortex glass phases has not been established 
conclusively. The boundaries shown between phases (I) and (II) and between phases (II) and (III) 
may actually correspond to rapid crossovers instead of true phase transitions. 
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FIG. 4. The diagonal correlation function g^^> {p, n = 0), obtained using the formalism described 
in section IIIC, is plotted as a function of the in-plane separation p/a, where a is defined by 
vra^ = ^q/Bq. The results shown are for disorder strength F' = 3 x 10"^F^ (see text for the 
definitions of F and F'), Bq = 0.3 T, and T = 14.5 K and 30 K. 
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FIG. 5. The off-diagonal correlation function g^'^'{p,n = 0), obtained using the formalism 
described in section IIIC, is plotted as a function of the in-plane separation p/a, where a is defined 
by TTO? = ^o/Bq. The parameter values are the same as those in Fig. 4. 
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FIG. 6. Fourier components pc of the periodic density field at the freezing temperature, as cal- 
culated in the density functional theory (triangles), plotted versus |Gp on a semi-logarithmic scale. 
The solid line corresponds to a Gaussian form for pc (see text) which has the same Debye- Waller 
factor as that obtained in the density functional theory at the freezing transition. 
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FIG. 7. Magnetic field distribution function n{B) in tlie crystalline phase at the freezing tem- 
perature for (a) Bq = 0.3 T and (b) Bq = 0.05T. The thermally broadened time-averaged density 
distribution in the crystalline state is obtained using the results of the density functional the- 
ory. The distribution n(B) for a slightly broadened Abrikosov flux-lattice (dashed lines) is also 
shown for comparison. Note that the distribution given by the density functional calculation is 
considerably narrower and more symmetric. 
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FIG. 8. The second moment [AB'^] of the magnetic field distribution n{B) versus temperature 
T for Bq = 0.3 T, obtained using the results of density functional theory for the Fourier compo- 
nents of the density field and a Debye- Waller approximation for the low-temperature regime. For 
comparison, the results for a perfect Abrikosov lattice are also shown (solid line). Corrections due 
to a finite core size (see text) are included. Note the curvature of the moment as a function of T, 
which is opposite to that obtained for a perfect Abrikosov lattice assuming a two-fluid form for the 
temperature dependence of A. 
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FIG. 9. The third moment [Ai?^] of the magnetic field distribution n{B) versus temperature 
T for Bo = 0.3 T obtained using the results of density functional theory for the Fourier com- 
ponents of the density field at the freezing transition and a Debye- Waller approximation for the 
low-temperature regime. Corrections due to a finite core size (see text) are included. 
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